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Abstract 

We investigate the analytic structure of the Borel-summed perturbative QCD amplitudes in 
the complex plane of the coupling constant. Using the method of inverse Mellin transform, we 
show that the prescription dependent Borel-Laplace integral can be cast, under some conditions, 
into the form of a dispersion relation in the a-plane. We also discuss some recent works relating 
resummation prescriptions, renormalons and nonperturbative effects, and show that a method 
proposed recently for obtaining QCD nonperturbative condensates from perturbation theory is 
based on special assumptions about the analytic structure in the coupling plane that are not valid 
in QCD. 
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I. INTRODUCTION 



The QCD amplitudes have a complicated analytic structure in the complex plane of the 
coupling constant, a s = g 2 . As proved in , the infinitely many multiparticle branch 
points at large energies result, via renormalization group invariance, in an accumulation 
of essential singularities near the origin a s = 0. Since in massless QCD the multiparti- 
cle hadronic states are generated only by a nonperturbative confinement mechanism, these 
singularities can show up only beyond perturbation theory. However, the perturbative am- 
plitudes themselves are expected to have a complicated structure of singularities, due to the 
fact that the perturbation series is divergent and Borel non-summable. The presence of the 
renormalons on the real axis of the Borel plane induces singularities of the amplitudes as 
functions of the coupling constant. 



In a recent paper [3] , arguments based on analyticity in the coupling complex plane were 
used to suggest the possibility of calculating genuine nonperturbative quantities, like QCD 
condensates, from pure perturbation theory. The analytic structure and its connection with 
infrared renormalons were further discussed in . Motivated by this recent interest in the 
problem, we investigate in the present work the analytic structure in the complex coupling 
plane of the Borel-summed amplitudes in perturbative QCD. We use the mathematical 
techniques applied in (f| Q, sjl , which allow us to express the Borel integral as a dispersion 
relation in the coupling plane. In the last section we shall make a few comments on the 
papers flfl jj. 

We consider for illustration the Adler function in massless QCD 



where Tl(Q 2 ) is the current-current correlation function calculated for euclidian arguments 
Q 2 > 0. It is known that the perturbation expansion of T>pt(cl) in powers of the renormalized 
coupling a = a s (Q 2 )/7c is divergent and not Borel summable (see [9| and references therein). 
The attempts of performing the summation by a formal Borel-Laplace integral encounter 
the difficulty that this integral is not well-defined. We consider the Borel transform B(u) 
defined in the standard way in terms of the perturbative coefficients d n of T>: 



B{u) = t d ^(^\ , (2) 
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FIG. 1: The Borel plane for the Adler function 

where /3q = (33 — 2rif)/l2 is the first QCD beta-function coefficient with rif the number 
of flavors. It is known, from the n\ large order growth of d n , that B(u) has singularities 
(ultraviolet and infrared renormalons) on the real axis of the it-plane j9|. For the Adler 
function, the ultraviolet renormalons are placed along the range u < — I and the infrared 
renormalons along u > 2 (see Fig. QJ. Due to the infrared renormalons, the usual Borel- 
Laplace integral is not well-defined and requires an integration prescription. Defining 

ooiie 

V ( pMa) = — e Po- B(u)du = — lim / e ^-B(u)du, (3) 

C± Q±ie 

one can adopt as prescription, for each value of a with Re a > 0, either T>pj.(a) or T>pj, (a), or 
a linear combination of them, with coefficients £ and 1 — £ such as to correctly reproduce the 
known low-order expansion of T>pt{cl) (which is obtained by truncating the Taylor expansion 
(J2J) at a finite order N). We consider in particular the principal value (PV) prescription 

V^Pia) = \[V^{a) + V#ta)\ . (4) 

Once a prescription is adopted, one has a well-defined function of a, different prescriptions 
yielding different functions. In the next section we shall study the analytic properties of 
these functions in the complex a-plane. 
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II. DISPERSION RELATIONS FOR THE BOREL SUMMED AMPLITUDE 



The analytic properties of the integrals Q with respect to the variable a can be studied 
with standard mathematical techniques. In the present work, we use the method of inverse 



Mel 
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in transform, applied for the first time in the context of Borel summation in QCD in 

— i n 

12j. For details of the mathematical procedure used below, we refer to [6], where the 
same method was applied for investigating the momentum-plane analyticity structure of the 
Adler function in the large-/?o limit. 

The inverse Mellin transform of the function B(u) is defined by [l^ 

uo+ioo 

w(t) = — [ B(n)r B - 1 d!i (5) 

uo—ioo 

and admits the inverse relation 

oo 

B{u) = [w(t)t- u (1t, (6) 







which gives B(u) in the strip — 1 < Rew < 2 parallel to the imaginary axis |l3j. The above 



relations are valid if the Borel transform satisfies the condition 



u +ioo 

— J \B(u)\ 2 du<w, (7) 

uq—ioo 

where uo is a point located on the real axis, between the branch points, — 1 < uq < 2 (see 
Fig.l). This condition strongly restricts the asymptotic behaviour of the Borel transform, 
and it is not known whether it is obeyed or not in QCD. The condition (J7J) is however 
satisfied in some particular cases of physical interest. One example is the summation of a 
chain of diagrams in the large-/5o (or large-ri/ limit), which leads to the Borel transform 



14. 



15|: 



BM=™± a , (8) 

3(2 -«) t 2 [P-(l-n) 2 ] 2 ' 

working in the ^-scheme, where all the exponential factors are included in the definition of 
the coupling. As a second example we take the case of a finite number of renormalons with 
branch-point singularities, in particular the contribution to B(u) of the leading infrared and 
ultraviolet renormalons 

K K' 

(n) ~ (i - u/p)-^ + (i + u/p'Y +1 ' 



where p > 0, p' > and the constants K and K' represent the strength of the corresponding 



singularities. We note that 



were calculated 



in 



and 



or the Adler function p — 2, p' — 1, and the exponents v and v' 



171 ]. respectively. 



The function w{j) defined in (j5J) can be calculated by closing the integration contour 
along a semi-circle at infinity in the u plane and applying the theorem of residues. For 
|t| < 1 the contribution of the semi-circle at infinity vanishes if the contour is closed in 
the right half of the u plane, while for |r| > 1 the contour must be closed in the left half 
plane. Therefore one obtains different expressions for the distribution function at |r| < 1 
and |t| > 1. By Cauchy's theorem these functions, which we denote by w < (t) and w>(r), 
have the representations 



■UMT 
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/ B(u) t u ~ du 



(10) 



1 

2~ni 



J B(u) r^du, 



the integration (Hankel) contours Cir and Cuv being indicated in Fig. 1. 



The inverse Mellin transform w was calculated in Ref. 
large-/5o limit, when the Borel transform has the expression 



ll| for the Adler function in the 
. In this case 



w > (t) 



lnr 



1 + rl 



L 2 ( 
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1 + rl 



Ln(-T- 1 ) -lnrlnfl + r' 1 ) 



where L 2 (x) = — Jq y ln(l — t) is the dilogarithm. The physical interpretation of w as the 
distribution of the internal gluon virtualities in Feynman diagrams was also pointed out in 

Q. 

The function id can be also calculated explicitly for a finite number of renormalons with 
branch-point singularities. For instance, for the Borel transform written in we obtain 
3: 



K 



W < (T 



w>{t) 



K' 



(pY +1 



lnr) 



7+i ( m r ) 



(12) 



T(v' + 1) " ' rP' 

We now proceed to the evaluation of the integrals 0, taking a in the right half plane, 
Re a > 0, where we assume that they converge. Our aim is to obtain a representation of 



Upr(a) in terms of the inverse Mellin transform w. To this end we rotate the integration 
contours C±, whitout crossing singularities, up to a line parallel to the imaginary axis, where 
the representation © is valid. It is easy to check that for a in the upper half plane (Im a > 0), 
the contribution of the quarter of the circle at infinity vanishes if the rotation is performed 
in the upper half of the w-plane. Using the representation (jHJ, valid along the imaginary 
axis, and performing the integral with respect to u (for details see [(|), we obtain 

P&(a) = T / 1 xi dT ' Ima >°- (13) 

Pd J ~a r m T 

Similarly, for a in the lower half plane (Ima < 0) the integration axis in the expression of 
Vp^(a) can be rotated in the lower half of the w-plane, up to the negative imaginary axis, 
leading to 



»»W = ft/ i utt7 dT ' Ima<0 ' (14) 

We notice further, recalling the definition (jSJ), that the first relation (fTUj) can be expressed 



as 



DjJ(a) = V$(a) + 2ia < (a) (15) 
where we introduced the notation 

ff (o) = ^[tw(t)] _ 1 . (16) 

In particular, using the expression (|12|) of w, we obtain for the leading renormalons: 

a<(a) = K -f^— £ e"^ (fto)-" , (17) 

^ - K ' W^vA^ ( ~ ° ar ' ■ 

It is important to emphasize that the relation (|TH|) is valid for Re a > (or equivalently for 
\t\ < 1), i.e. in the whole right-half of the a-plane. 

The relations (|13|). (|14|) and (fl3j) are the basis the derivation of the dispersion relations 
given below. We first note that, unlike the original representations © which converge 
only for Re a > 0, the representations |T3|) and (|T4|) can be analytically continued in the 
corresponding upper (lower) half of the complex a plane, outside the real axis. Moreover, 
it is easy to convert them into a dispersion representation in the variable a. We first split 
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the integral in two integrals, one from to 1 (where id = w<), and the other from 1 to oo 
(where w = u>>), and perform in each interval the change of variable 

1 da' 

lnr = - ^' ~ = M^- (18) 

Using the relations (fT3j) - (jl4|) thus transformed, together with ()15|) and (j!6|) . we finally express 
T>$(a) as 

^f+l/ x a, f er<(a') da' a r a > (a')da' _ 

Pff a = -/ +- / v Ima > ' 19 

7r J a (a' — a) 7r j a'(a — a) 

v 7 -oo v ' 

^r+i/ \ a / > cr < (a')da' a /• cr>(a')da' , . 

pW a = - / <V T + - / \ + 2»<7<(a , lma<0. 

7r j a a — a) 7r j a (a — a) 

From the definition and the properties of the inverse Mellin transform w, it follows that 
the function cr<(a) can be analytically continued in the complex a-plane. This property is 
seen explicitly in the case of one infrared renormalon in (|17j) . Thus, the expressions (J19|) are 
analytic functions in the upper (lower) half of the complex a-plane, outside the real axis. 
From (j!7|) it is seen that cr>(a) is real for a < 0, while cr<(a) is real for a > 0. Therefore, 
the spectral functions of the dispersion integrals are real. 

Taken together, the dispersion relations (fTTIJ) define a single analytic function, Dp^(a), in 
the whole cut a-plane. This function may have a discontinuity across the positive axis, due 
the Cauchy dispersion integral along a > and the additional term 2ia < (a) in the second 
relation. A simple calculation shows however that 

, r ov(a')da' a °r cr > (a')da / „ . . 

lira pW fl ± ie ) = -P / <\ > + - / >\ > + i a< (a), a > , (20) 
e->o+ 7r J a a — a n J a (a — a) 

V ' -oo v ' 

where P denotes the Cauchy principal value. This relation shows that the function "Dp^(a) 
is well defined along the positive axis, but has there an unphysical imaginary part equal to 

ImP^(a) = t7<(a), a>0. (21) 

From (fE?j) it follows that T>^ (a) has actually a discontinuity along the negative axis, given 
by 

pfj (a + ie) - V$ (a - ie) = 2i [o> (a) - a< (a - ie)} , a < , (22) 

in terms of the real spectral function cr> and the analytic continuation of <r< up to the lower 
edge of the negative semiaxis (where, as seen from (jl7j) . it is complex). 



The function Vpp(a) satisfies a dispersion relation similar to ([19)1 . with an additional 
term — 2ia < (a) in the right hand side, as follows from (fTB^) . As above, one can show that 
^pr( a ) * s well-defined for a > 0, but it assumes there complex values. The unphysical 
imaginary part is eliminated if we adopt the principal value prescription defined in (0J), for 
which we obtain: 

^(pv), s a f cr < (a') da' a °r a > (d)dd , . . . 

V " W = + - / d^V— <(«). M«>0 (23) 



-oo 

7r j a (a — a) n J a [a — a) 



o 

For a on the positive semiaxis, is easy to check that 

(PV) f ■ \ a n f cr < (d)dd | a °f a > (a')da' 



lim V^\a + ie) = -P / <\ ' . + - / f// . , a > . (24) 
e^o+ ^ J ' ix J did - a) n J did - a) 

o v ' -oo v ; 

Therefore, T>^p^\d) is well-defined and real on the positive real semiaxis, as required by 
general principles. The expressions (|23|) define an analytic function which satisfies the reality 
condition Vp P ^\d) = [£>^ / ' ) (a)]* in the whole complex a-plane cut along the negative 
semiaxis, where it has a discontinuity 

V ( pp (a + ie) - V { pp (a - ie) = 2i [a> (a) - Re a > (a)] , a < 0. (25) 



As argued in [6|, this discontinuity vanishes only under strong restrictions on the asymptotic 
behavior of B(u), which are satisfied neither in the simple cases considered here, nor, most 
probably, in full QCD: namely, the L 2 condition (JJJ) must be satisfied not only by B(u), but 
also by the product B(u) shi7ru (for technical details, see Ref. [(|). 

We point out that the imaginary part of Dp^ (a) along the positive semiaxis vanished 
due to a precise cancellation of the imaginary part of the integrals and the last terms in Eqs. 
f!23|) . It is easy to check that for a general linear combination of Vp^, with coefficients £ 
and 1 — £ as discussed above (0J), this cancellation no longer holds. Therefore, the principal 
value prescription is the most suitable choice if one wants to preserve in perturbative QCD 
the analytic properties of the true amplitudes. 

In this section, we obtained the nontrivial result that the inverse Mellin transform can be 
used to derive from (jHJ) dispersion relations for the perturbative amplitudes V P ^ r )(a) J which 
allow an analytic continuation of the Borel integral into the left-hand half-plane Re a < 
and explicitly exhibit the singularities and the discontinuities across the cuts. 



The results obtained here, besides expressing the Borel integral (j3J) in the more suitable 
form of a dispersion relation, will be useful in discussing the validity of some assumptions 
on analyticity in the a plane made in the literature, as we show in the next section. 



III. COMMENTS 

In this section we shall make a few comments on the recent papers 0, |4|, |^|, related to 
the present work. 

In Ref. ji] the author considers, in the case of a single infrared renormalon, the problem of 
removing the unphysical imaginary part of T)p^ (a) by subtracting a suitable regularization 
function from it. Our results in this particular case are consistent with Indeed, Eqs. 
(fT9]l and (|2*H|) show that the function V^p^ (a) is obtained by subtracting from T>pj, (a) the 
function 

A(o) = z«7<(o) , (26) 

which, using ((17)). can be written for a = \a\e 1 ^ as 

2^+i / n \ 2_ 

A(a) = K— — -He "°« \(3 a\- u [sm tfju + i cos ipu] . (27) 

This expression coincides with Eq. (17) of Q, derived using arguments based on regularity 
with respect to the parameter v. 

The analytic structure in the coupling constant plane was recently considered also in Ref. 
0| , where arguments based on analyticity were used in support of the claim that the QCD 
condensates can be determined using the coefficients of the perturbation series. In what 
follows we shall briefly analyse the validity of this claim. 

The author of jjj uses the analogy with some semiclassical models [19], where a specific 
contribution (for instance, multi-instantons) is regular for negative couplings and can be 
obtained for positive couplings by analytic continuation. He invokes the heuristic argument 
according to which a perturbation series with a sign-nonalternating, n\ large order behavior, 
can be summed at a < by a Borel integral along the negative axis in the -u-plane, where 
it becomes sign alternating. We note however that sign alternation does not necessarily 
imply Borel summability: this property is violated if a nonalternating component is present, 
however negligible it may be in comparison with a strong sign-alternating component of the 
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series. This is exactly the situation in QCD: due to the ultraviolet renormalons, the Borel 
integral along the negative u-axis is not well defined. 

Another conjecture adopted in [3] is related to the nonperturbative amplitude V NP (a), 
which must be added to T>px{a) in order to compensate its unphysical imaginary part: 

lmV NP (a + ie) + lmV PT (a + ie) = 0, a>0. (28) 

According to current interpretations , the cancellation is expected to occur if both terms in 
the above relation are calculated with the same prescription. In |3[ the author supplements 
PS]) by a specific assumption about the nonperturbative amplitude, taking it of the form 

V NP {a) = Ce-W a (-f3 a)-\ (29) 

where v is the branch-point exponent of the infrared renormalon in Q and C is related to 

n 

the gluon condensate (we use the parametrization given in |5[, with our notation a = a s /n). 
This expression is regular for negative a, having a discontinuity only on the positive semiaxis 
in the a-plane, where 

2_ 

V NP (a±ie) = C e <V ((3 a)~ u [cosnu =p isimru] , a>0. (30) 

Using this relation and the imaginary part of the perturbative amplitude given in (}2~Tj) and 
(JT7J), condition (|2~8*j) gives the relation 

K = C^-^- } (31) 
2^+i T{-v) V ; 

from which, according to Q], one could obtain the nonperturbative parameter C (the gluon 
condensate), using the strength K of the infrared renormalon computed from the perturba- 
tion series. 

The relation (}3"T]) implies that K vanishes for nonnegative integer u, when the renormalon 



1/(1 — u/2) 1+u becomes a pole. It is known however that poles are actual 



y obtained from 



some chains of Feynman diagrams, in the large-/3o (or large-ri/) limit Q, [l3- In H the 

author discusses this limit, taking the exponent v of the form v = k + x/ '/3q + ■ ■ ■, with k an 

integer. In this limit, the factor T(—v) ~ flafx m ^ e denominator of (J31)) is compensated 

by the factor (3q in the numerator, leading to a finite nonzero limit for K, if C tends to a 

nonzero constant in the large-/?o limit. But then Eq. (|3"Uj) implies that the real part of the 

rn 

nonperturbative amplitude T>np is nonvanishing in this limit (we keep, following [5|, the 



product fioa constant, which is legitimate, as seen in particular in the one-loop expression 
a = p olri Q2/ A 2 )- On the other hand, the perturbative amplitude V PT is subleading for large 
/3q, due to the factor 1//3q in front of the Laplace-Borel integral (JHJ). A choice of the constant 
C of the same order, i.e. C ~ c/ /%, would imply from ([HI]) that K vanishes when (3q — > oo. 
So, it follows from 0,0] that in the large- (3q limit either the total Adler function is dominated 
by the real part of the nonperturbative amplitude, or the renormalon residue vanishes. 

The above implications depend however, in a crucial way, on the conjecture made in [^J 
that the only singularities of the nonperturbative amplitude are along the positive semiaxis. 
Actually, singularities along the negative semiaxis (produced by the ultraviolet renormalons 
in the OPE coefficients) cannot be excluded. In |3j it is claimed that the influence of the 
ultraviolet renormalons can be suppressed by an appropriate conformal mapping jlfij . A 
numerical suppression, however, does not imply removal of the corresponding singularities. 
As is shown in joL even if B{u) is expanded in powers of the optimal conformal mapping of 
the Borel plane [20], the function Dprifl) does have singularities along the negative semiaxis. 
Therefore, the general expression of T>Np(a) may have singularities both for a > and for 
a < 0. To give an example, we add to (jHOJ) the analogous branch-point term suggested by 
Eq.flHl), taking 

V NP {a) = Ce'^ {-fa a)- u + C'e~ ^(/3 a)~\ (32) 

where C and C depend in general on v and C = C' R + iC'j is complex (as the condition 
(I28|) holds in a definite prescription, we specifically refer here to the prescription leading to 
Vpj,). Then condition (|28j) gives 

-Camnv + C'j + K /\ ^ =0, (33) 

1 [y + l) p 

which does not necessarily imply that K vanishes when v is a nonnegative integer. Moreover, 
the real part of the nonperturbative amplitude (|3*2*|) contains the additional parameter C' R 
not specified by the condition ()33|) . which, combined with the first term in the expression 
(|3*2*|) can make KeV^p subleading in the large- (3q limit. This counterexample disproves 
the conclusions (that appear to be inherent in the method of jjj O]) that the renormalon 
residue vanishes when the singularity is a pole, or the amplitude is dominated in the large-/9o 
limit by the nonperturbative part. Even more important, however, is the fact that the real 
part of the nonperturbative amplitude cannot be determined from the requirement (|28J) and 
analyticity arguments. The possibility, advocated in jjjj], of obtaining the gluon condensate 
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from the perturbation series fails to work here. 

In the present paper, by using the inverse Mellin transform, we showed that the Borel- 
summed QCD perturbative amplitudes J3J) satisfy dispersion relations which explicitly ex- 
hibit the singularities in the complex a-plane and the discontinuities across the cuts, and also 
allow the analytic continuation of the Borel integral into the left-hand half-plane Re a < 0. 
We saw that the derivation relies on some special hypotheses about the properties of the 
Borel transform in perturbative QCD: we namely assumed that the Borel transform has 
no singularities in the complex w-plane except for branch-points on the real axis, with a 
holomorphy gap around the origin, its asymptotic behavior being such that the inequality 
(J7J) holds and the inverse Mellin transform (jSJ) exists. We point out that the holomorphy 



of B(u) in the doub 
field theories 
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e-cut it-plane is expected on general grounds to hold in renormalisable 
. The validity of this condition in QCD is presently almost universally 
adopted as a plausible assumption. As for the asymptotic condition, its validity in full 
QCD is not guaranteed, but it is satisfied in the large-/3o limit and for a finite number of 
renormalons with branch-point singularities. 

In the second part of the paper we made some comments on the papers 0, U, 0], where 
arguments based on analyticity were used to discuss the connection between renormalons 
and nonperturbative quantities in QCD. We investigated some specific conjectures made in 
and showed that they are in conflict with the analyticity properties related to the 
ultraviolet and infrared renormalons, which throws doubts on the method of calculating 
QCD condensates from perturbative expansions. 
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